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F^ ' Abstract 

O I We consider the fractal dimensions of critical clusters occurring in configurations of 

a q-state Potts model coupled to the planar random graphs of the dynamical trian- 
gulations formulation of Euclidean quantum gravity in two dimensions. For regular 
lattices, it is well-established that at criticality the properties of Fortuin-Kasteleyn 
clusters are directly related to the conventional critical exponents, whereas the cor- 
^^ - responding properties of the geometric clusters of like spins are not. Recently it has 

,—1- I been observed that the latter are related to the critical properties of a tricritical 

^^ ] Potts model with the same central charge. We apply the KPZ formalism to develop 

a related prediction for the case of Potts models coupled to quantum gravity and 
employ numerical simulation methods to confirm it for the Ising case q = 2. 
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. >r , 1 Introduction 



The relation between the percolation problem and tliermal phase transitions 
of lattice spin systems has been a question of intense research for at least three 
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decades. Clusters of even spins are natural objects occurring in the analysis of 
phase ordering processes and nucleation [l| , and a theory of critical phenom- 
ena in terms of purely geometrical objects appears appealing. In this context, 
it had long been surmised that a continuous phase transition of a spin sys- 
tem might be accompanied (or, in fact, caused) by a percolation transition 
of the clusters of like spins {geometric clusters), the appearance of a per- 
colating cluster sustaining the onset of a non-zero magnetisation [2]. While 
for the special case of the Potts model in two dimensions it turned out that, 
indeed, the thermal phase transition point coincides with the percolation tran- 
sition of the spin clusters, this behaviour is not generic and does not occur in 
three-dimensional systems J3[. Also, the critical exponents associated to the 
percolation of geometric clusters are not directly related to the thermal expo- 
nents of the spin model [4]. However, a close relation between the percolation 
and thermal phase transitions can be established by considering stochastically 
defined clusters (or droplets) as they occur in the Fortuin-Kasteleyn (FK) rep- 
resentation of the Potts model, and it can be shown that, in fact, the Potts 
model is equivalent to a site-bond correlated percolation problem p such that 
the corresponding critical exponents agree. This identification of the proper 
cluster objects [FK clusters) percolating at the thermal phase transition sub- 
sequently also allowed for the design of cluster algorithms for the efficient 
simulation of Potts models in the vicinity of the ordering transition p, |7|, 
beating the observed critical slowing down of local update algorithms. Sim- 
ilarly, relations of continuous-spin models to percolation problems could be 
established and corresponding cluster algorithms formulated [7|, |8|, such that 
the continuous phase transitions of many standard models of statistical me- 
chanics are by now understood in terms of the percolation properties of some 
suitably defined ensemble of stochastic clusters. 

Although not in the universality class of the thermal phase transition, the 
clusters of like spins or geometric clusters still undergo a percolation tran- 
sition in the course of thermal phase ordering. This transition is in general 
not equivalent to ordinary (site or bond) percolation jJ] and it remains an 
interesting open question to determine the general critical behaviour of clus- 
ters of aligned spins. For the case of the two-dimensional Ising model, it has 
been conjectured and numerically verified that the geometric clusters are de- 
scribed by the g = 1 tricritical Potts model J9|; this correspondence can be 
understood from a direct construction starting from the dilute Potts model 
(lO|. Subsequently, analogous conjectures for the 2 < g < 4 Potts models were 
made and some of them substantiated by numerical simulations [ll|. Analyti- 
cal calculations concerning clusters occurring in systems of statistical physics 
and their boundaries, traditionally based on methods of conformal field theory 
and the Coulomb gas [l2], have recently seen major advances from the insight 
that fractal random curves can be described in a framework dubbed stochastic 
Loewner evolution (SLE) |l3|, |l^. Collecting these observations, a more sys- 



tematic analysis of the relation between the critical and tricritical branches of 



the Potts model and their connection to the FK and geometric clusters has 
been performed, resulting in the identification of exact values for the different 
cluster fractal dimensions and their numerical verification |l5lllq. 

Coupling spin models to the planar random lattices of the dynamical triangu- 
lations model of two-dimensional Euclidean quantum gravity |l7[ corresponds 
to the introduction of a particular type of annealed connectivity disorder. The 
resulting randomness is strong, leading to a change in critical behaviour of 
virtually all types of coupled matter variables |l8[. Geometrically, it is char- 
acterised by a large fractal dimension d^ ^ 4 oi these lattices of topological 
dimension two. It is interesting to see how the relation between geometric and 
FK clusters of the Potts model works out far away from the regularity of a 
Bravais lattice. The mentioned relations between geometric and FK clusters 
have not yet been rigorously established, such that evidence from further mod- 
els is highly welcome support. Besides, the behaviour of these fractal cluster 
objects on lattices which are themselves highly fractal is of particular interest 
in itself. 



2 Dynamical triangulations and the KPZ formula 



The dynamical triangulations approach provides a constructive model for Eu- 
clidean quantum gravity in general dimensions [l7[. It regularizes the path 
integral over fiuctuating metrics naturally occurring in an attempt to quantize 
the gravitational interaction by a sum over combinatorial manifolds realised as 
simplicial complexes [19| . In two dimensions, the resulting canonical ensemble 
of discrete surfaces can be defined as that of all possible gluings of a given 
number N2 of equilateral triangles to a closed surface of fixed (usually spher- 
ical) topology, where all resulting triangulations are counted with the same 
weight in the partition sum. This is a purely combinatorial definition, and 
due to the equilaterality of the triangles the resulting triangulations cannot in 
general be embedded in the Euclidean M^ and, in fact, one is only interested 
in their intrinsic geometry. Many counting problems related to these graphs 
can be solved exactly by using matrix integrals or general combinatorial tech- 
niques J20|. For our purposes it is sufficient to note that the resulting random 
graphs are highly fractal with blobs ( "baby universes" ) of arbitrary size be- 
ing connected to the main graph body with a minimal number of links. This 
structure entails an internal Hausdorff dimension dh = 4 |l7| . Decorating the 
vertices (or edges, faces) of the graphs with matter variables, corresponding 
to an annealed geometrical average, leads to a change of universality class at 
criticality, expressed through a dressing of the conformal weights A of the 



matter fields given by the KPZ formula [1 

VI - c + 24A - v^r^c 
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where c denotes the central charge of the coupled matter model. Analogously, 
the change of the string-susceptibility exponent 7s can be expressed in terms 
of c ET 



As mentioned above, the fractal properties of FK and geometric clusters of 
the square-lattice Ising model have been studied rather extensively |9LllCill6l|. 
Here, we are interested in the (normalised) fractal dimensions d^ /d resp. 
d^/d of the incipient percolating Fortuin-Kasteleyn resp. geometric cluster at 
criticality (where d = 2 denotes the spatial dimension). For the FK clusters. 



d^/d = 1 — P/du = 15/16 is an exact result 2j|. From the identification 
of the geometric Ising clusters with the FK clusters of a g = 1 tricritical 
Potts model, one finds d^/d = 187/192 [lO|. These values are related to the 
conformal weights A above as Ac = 1 — dc/d 2^. To describe the coupling 
of the two-dimensional Ising model to quantum gravity, the resulting weights 
A^^ = 1/16 and A^ = 5/192 have to be dressed according to Eq. (1). Since 
for both, the Ising and g = 1 tricritical Potts models, the central charge 
c = 1/2, from (1) we find A^^ = 1/6 and A^ = 1/12. Therefore, we expect the 
following cluster fractal dimensions for the Ising model coupled to dynamical 
triangulations, 

rf^^/4 = 1 - A^^ = 5/6, Jg/4 = 1 - Ag = 11/12, (2) 



which now have been written in units of the Hausdorff dimension dh of the 
graphs coupled to the spin model, the value of which is not known exactly but 
numerically found to be consistent with dh = ^ for all < c < 1 23|, |24 . 



3 Numerical simulation method and results for the Ising model 



For a numerical simulation of the system, due to the annealed nature of the 
disorder both, the underlying dynamical triangulations as well as the coupled 
Ising spins, must be updated in parallel. An ergodic set of Monte Carlo updates 
for the dynamical triangulations is given by the so-called Pachner moves J25|. 
For the canonical ensemble of a fixed number N2 of triangles in two dimensions, 
these reduce to the following flip between two adjacent triangles. 





where the dashed hnes indicate the effect on the dual cfr" graph. To improve 
the decorrelation of adjacent configurations in the Monte-Carlo Markov chain, 
this local update is complemented by intermittent non-local rewirings of the 
blob structure of the graph known as "baby-universe surgery method" [2&..27|. 
We consider the case of non-degenerate triangulations, corresponding to dual 
one-point irreducible 0^ Feynman diagrams and place the Ising spins on the 
faces of the triangulation or, equivalently, the vertices of the dual 0^ graphs. 
The coupled Ising model is being updated according to the Swendsen-Wang 
cluster algorithm [(| in order to alleviate the critical slowing down of the 
spin variables expected since simulations are performed at the exactly known 
inverse critical temperature /3c = | In ^ :^ 0.7733 of the coupled system J28|. 

Although the fractal dimensions d^ and d^ could be determined directly by 
means of an appropriate geometrical analysis of the clusters, it is convenient 
to exploit relations to more easily accessible quantities. The fractal dimension 
of Fortuin-Kasteleyn clusters, d^, is related to the magnetic susceptibility 
exponent 7 = 7^^ as dc/dh = |(1 + 'y/dhi') J29|, where u denotes the critical 
exponent of the correlation length. The susceptibility y = x^^ can be sampled 
by the following estimator in the cluster language J30| , 




x-TT E-. =ir Eic^-n. (3) 



where cxj = ±1 denotes the Ising spins. The sum in the second expression is 
over all FK clusters Cj (including the percolating cluster as well as isolated 
sites) of a given configuration and \Ci\ denotes the number of spins in cluster i. 
Close to the critical point and for a system of finite size A^2; standard finite-size 
scaling (FSS) arguments suggest that x scales as 



Thus, Eqs. (3) and (4) allow for a FSS determination of the fractal dimension 
dc^/dh of the FK clusters from information about the cluster distribution 
naturally produced by the Swendsen-Wang cluster update of the spins. In 
complete analogy, a "geometric susceptibility" x^ can be estimated with Ci 
now denoting the geometric clusters in relation (3), and its FSS is described 
by Eq. (4) with exponent d^/dh. 
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Fig. 1. Finite-size scaling of the magnetic ("FK clusters") and geometric ("ge- 
ometric clusters") susceptibilities x of the Ising model coupled to regular dy- 
namical triangulations with A'^2 triangles and at the critical inverse temperature 
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0.7733. The lines show fits of the functional form (4) to the data. 



As a check and gauge for the method, we determined the cluster fractal di- 
mensions for the Ising model on the square lattice, using L x L lattices with 
L = 8, . . . ,512. From fits of the form (4) to the (magnetic and geometric) 
susceptibility data, we find d^^/d = 0.93 756(15) and d^/d = 0.97370(13) 
(including lattices of sizes L = 16, . . . , 512), in very good agreement with the 
exact results dl^/d = 15/16 ^ 0.93 750 and d^/d = 187/192 ^ 0.97396. 
For the dynamical triangulations model, only graphs of somewhat smaller 
sizes can be equilibrated, and we consider volumes N2 = 256, 512, . . . , 65 536. 
Figure 1 shows the resulting susceptibility data in a doubly logarithmic plot 
together with fits of the form (4) to the data. Due to rather strong finite-size 
corrections, only the largest lattice sizes can be included in the uncorrected 
fit (4) with reasonable fit quality. For the range N2 = 8192, . . . , 65536, we 
arrive at the estimates dl^/dh = 0.84 301(93) and rfg/4 = 0.92 220(55), in 
marginal agreement with the conjectured values of Eq. (2), d^/dh ~ 0.83 333 
and d^/dh ~ 0.91 667. From previous experience with the dynamical trian- 
gulations model, one indeed expects very strong finite-size corrections to be 
present J2J|. These result from the small effective linear extents of lattices of 
a given number of triangles reflected in the large fractal dimension dh ~ 4. As 
dominant contribution, we expect analytic corrections to the effective linear 
extent, 



L,s{N2 



LqNI''''^ + Li + L^N^^''"'^ + 



(5) 



such that the susceptibihty should scale as x ~ [Lcq{,N2)\^^'^ . Since the data 
are not precise enough to determine an additional exponent, we fix d/i = 4 in 
the expression for Les{N2). Already with one correction term only (L2 = 0), 
the fit results, 



rf^^/4 = 0.8291(19), c/g/4 = 0.9132(12) (6) 



for the range N2 = 1024, . . . , 65536, move considerably closer to the con- 
jectured values. With both correction terms (i.e., variable L2), we get full 
consistency with estimates J™/4 = 0.8349(118) and rfg/4 = 0.9164(68). 
Conversely, if we fix d^/dh resp. d^/dh at the conjectured values, the fits 
with one correction term give reasonable fit qualities with Q = 0.20 resp. 
Q = 0.03, and very high-quality fits are reached when including both correc- 
tion terms with Q = 0.66 resp. Q = 0.42. 



4 General results for the Potts model 



To understand the general relation between FK and geometric clusters and 
the critical and tricritical branches of the Potts model coupled to dynamical 
triangulations, we consider the following parametrisation of the critical g-state 
Potts model 0, 

a/^ = — 2cos(7r//i;), k, = , ?ti = 1, 2, 3, . . . , (7) 

m 



where 1 < k < 2 and the central charge c = 1 — 6/171(771 + 1), such that k = 2, 
3/2, 4/3, 6/5, 1 corresponds to the g = 0, 1, 2, 3, 4 Potts models with c = —2, 
0, 1/2, 4/5, 1, respectively. The conformal weights of the primary operators 
follow from the Kac table 22l |. 



\(m + l)r — ms]'^ — 1 
A = .^-^ ■ - - rs = 12 777 — 1 

4777(777 + I) 



and they are identified with the physical operators of the Potts model as 
follows, 



Ae = A2,l 



_3k 1 

~ T~ 2' 

A A 1 Sk 1 
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(9) 
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where A^ and A^- = Ac correspond to the leading energetic and magnetic 
operators, respectively, Ah denotes the weight corresponding to the cluster 
hull, Aep the weight of the external perimeter and Arb the weight of the 
"red bonds" of the cluster, cf . Ref . [la| . The various fractal dimensions of FK 
clusters are given by the corresponding renormalization-group eigenvalues, 
da = Va = dh{l — Aq), where a = C, H, EP, RB. 

In Ref. jla] it was argued that in general the fractal dimensions of the geomet- 
ric clusters of the critical Potts model are identical to the fractal dimensions 
of the tricritical Potts model of the same central charge, to be reached by per- 
forming the central-charge conserving "duality transformation" k —>■ 1/k in 
Eq. (7). For the tricritical model, the leading conformal weights are identified 
as 

A, 
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(10) 
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which, as can be seen, directly follow from those of Eq. (9) by replacing k 

1/k. 



For the weights of the unitary minimal models of Eq. (8), the KPZ formula 

(1) can be written as 

Ar,s = ^il-K+\S-Kr\), (11) 



explicitly revealing that all weights of the minimal models dressed for the 
coupling to quantum gravity are rational numbers. Dressing the weights (9) 
of the critical branch yields 
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whereas the weights (10) of the tricritical branch become 

A. 



3 


'-^-\- 


K 


3 K 

4~ 2' 


Ah = -, 





Ac = 7- 2' ^H = -, (13) 

Aep = 1--, Arb = - + -. 



Note that now, in contrast to the regular lattice case of Eqs. (9) and (10), 
all weights depend linearly on the parameter k. The resulting FK and geo- 
metric cluster fractal dimensions for the Potts model coupled to dynamical 
triangulations are summarised in Table 1. 



5 Conclusions 



We have shown how the fractal dimensions associated with the geometric 
clusters of the Potts model coupled to the dynamical triangulations model of 
Euclidean quantum gravity in two dimensions can be inferred from a mapping 
to the tricritical branch of the Potts model. To allow for an application of the 
KPZ framework, the conformal weights associated to the fractal dimensions 
dc of the cluster, d^ of the cluster hull, ciEP of the external perimeter and c^rb 
of the cluster red bonds have been identified from the Kac table. Lifting the 
corresponding weights of the critical and tricritical branches to the dynamical 
triangulations results in rational values for the fractal dimensions of FK and 



Table 1 

Cluster fractal dimensions of the g-state Potts model on dynamical triangulations. 
The values of d-Ep/dh and d^^/dh for the geometric clusters are only formal since, 
by definition, dEp/d^ < du/dh- Thus one has d^p/d^ = d^/d^ and the geometric 
clusters are multiply connected. 
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geometric clusters for the g = 0, 1, 2, 3 and 4 state Potts models summarised 
in Table 1. The duality symmetry k, ^ 1/k present between the critical and 
tricritical weights on regular lattices, Eqs. (9) and (10), is no longer present 
in the dressed weights of Eqs. (12) and (13). As a peculiarity, we note that 
the fractal dimension of the geometric cluster hulls is a constant d^/dh = 1/2 
independent of q. A scaling or duality relation between the dimensions of the 
external perimeter and cluster hulls found valid for the regular case 



{dEp/d--){du/d- -) = — , 



(14) 



does not apply to the model coupled to quantum gravity, but instead one gets, 
(dEp/dh - \){dn/dh - ^) = 0. (15) 

Also, while the four fractal dimensions fulfil the identity dc—d^i = |((iEP— c^rb) 
for regular lattices [16| , for the random-lattice case we instead find dc — du = 
|(c^EP ~ c^rb)- a couple of further similar scaling relations can be formulated. 

Our simulation results for the case of the q = 2 Potts or Ising model cou- 
pled to dynamical triangulations show full consistency of the measured fractal 
dimensions of the Fortuin-Kasteleyn and geometric clusters with the predic- 
tions resulting from the KPZ mapping, confirming that the properties of the 
geometric Potts clusters on dynamical triangulations are described by the cor- 
responding tricritical Potts model of the same central charge. Corrections to 



10 



scaling are found to be much stronger for the random-graph model than for 
the square-lattice simulations performed as a gauge. These corrections, known 
to result from the small effective linear extents of the graphs due to their large 



Hausdorff dimension d^ ^ A [2J|, have to be explicitly taken into account to 
find consistency with the scaling predictions and satisfactory quality of the 
fits. 

It would be interesting to see whether, as expected, the fractal dimensions 
of the geometric clusters of the q ^ 2 Potts models coupled to dynamical 
triangulations follow the predictions summarised in Table 1, in particular for 
the case of g = 4, where the critical and tricritical branches coalesce and the 
KPZ mapping (1) becomes marginal due to central charge c = 1. 
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